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ABSTRACT 

Pre-protostellar cores likely represent the incipient stages of low-mass (~ IMq) star formation. Lynds 
1498 is a pre-protostellar core (PPC) and was one of the initial objects toward which molecular depletion 
and differentiation was detected. Despite the considerable scrutiny of L1498, there has not been a 
extensive study of the density and temperature structure as derived from radiative transfer modeling 
of dust continuum observations. We present deep SCUBA observations of L1498 at 850 and 450 /im, 
high resolution BEARS maps of the N 2 H + 1 — > transition, CSO observations of the N 2 H + 3^2 
transition, and GBT observations of the C3S 4 — > 3 transition. We also present a comparison of derived 
properties between L1498 and nearby PPCs that have been observed at far-infrared and submillimeter 
wavelengths. The L1498 continuum emission is modeled using a one-dimensional radiative transfer code 
that self-consistently calculates the temperature distribution and calculates the SED and intensity profiles 
at 850 and 450 /im. We present a more realistic treatment of PPC heating which varies the strength of 
the ISRF, Si sr f , and includes attenuation of the ISRF due to dust grains at the outer radius of the core, 
Ay- The best-fitted model consists of a Bonner-Ebert sphere with a central density of 1 — 3 x 10 4 cm~ 3 , 
R 0.29 pc, 0.5 < Sisrf < 1, A v w 1 mag, and a nearly isothermal temperature profile of « 10.5 K for 
OH8 opacities. C3S emission shows a central depletion hole while N 2 H + emission is centrally peaked. 
We derive a mean N 2 H + abundance of 4.0 x 10 -10 relative to H 2 that is consistent with chemical models 
for a dynamically young yet chemically evolved source. The observed depletions of C3S and H 2 CO, the 
modest N 2 H + abundance, and a central density that is an order of magnitude lower than other modeled 
PPCs suggests that L1498 may be a forming PPC. Our derived temperature and density profile will 
improve modeling of molecular line observations that will explicate the core's kinematical and chemical 
state. 



1. INTRODUCTION 

Pre-protostellar Cores (PPCs or starless cores) are 
thought to represent the initial stages of low-mass (M < 
few M ) star formation. PPCs are identified with dense 
molecular cores (n > 10 4 cm -3 ) which lack evidence for an 



internal heating source. Their column density is traced by 
dust continuum emission (Ward-Thompson et al. 1994; 
Ward-Thompson, Motte, & Andre 1999; Shirley et al. 
2000; Tafalla et al. 2002) and nitrogen-bearing molecular 
tracers such as NH3 and N 2 H + (Myers & Benson 1983, Ji- 



1 Jansky Postdoctoral Fellow at the National Radio Astronomy Observatory. The National Radio Astronomy Observatory is a facility of the 
National Science Foundation operated under a cooperative agreement by Associated Universities, Inc. 

2 2003 REU Summer Student at NRAO Socorro. 

3 2004 REU Summer Student at NRAO Socorro. 



2 



jina et al. 1999, Caselli et al. 2002, Crapsi et al. 2005). It 
is important to characterize their physical structure (den- 
sity, temperature, kinematical, and chemical) to constrain 
the initial conditions of low-mass star formation. It is 
also important to understand the evolution of the physi- 
cal structure of PPCs to trace the evolution of a nascent 
protostar. 

Lynds 1498 is a nearby molecular core (D — 140 pc) 
located in the west-central portion of the Taurus molec- 
ular cloud (see Cambresy 1999). L1498 was classified as 
a Pre-protostcllar core by its submillimctcr dust emission 
and lack of an IRAS detection at 100 /im (Ward Thomp- 
son ct al. 1994). The relatively weak submillimeter dust 
emission indicated that L1498 was a low density core. 

L1498 has been extensively studied with molecular 
line observations and using theoretical chemical modeling 
(Rawlings et al 1992, Taylor et al 1996). L1498 was ini- 
tially identified kinematically as a core potentially on the 
verge of protostellar collapse (Zhou et al. 1994). Sub- 
sequent molecular line observations have shown dramatic 
evidence for the freezing of gas phase molecules onto dust 
grains (e.g., H 2 CO, CS, CCS, C 18 0; Wang 1994; Lemme 
et al. 1995; Kuiper, Langer, & Velusamy 1996; Willacy, 
Langer, & Velusamy 1998, Tafalla et al. 2004); as a result, 
many molecular species cannot reliably trace the densities 
in the center of the core. Prominent exceptions are the 
nitrogen-bearing molecules NH 3 and N 2 H + which resist 
depletion (Flower, Pineau des Forets, & Walmsley 2005). 
Even though NH 3 and N 2 H+ remain in the gas phase in 
cold, dense molecular cores, the optical depth, gas temper- 
ature, and abundance profile must be reliably determined 
in order to accurately calculate the density structure of 
the core. 

Optically thin dust emission at submillimeter wave- 
lengths is a good tracer of temperature and density since 
the specific intensity is proportional to the integral of den- 
sity, opacity, and Planck function along a line-of-sight 

2^mnhv 3 f°° K u (s)n(s) 
I v = 5 / t r as (1) 

C 7 ° «P 

(Adams 1991, Shirley et al. 2003). Knowledge of the tem- 
perature and density structure is essential to a correct in- 
terpretation of chemical and kinematic models of L1498. 
Previous dust continuum studies of PPCs used radiative 
transfer models to constrain the density and temperature 
structure of the cores (e.g. Evans et al. 2001, Zucconi et 
al. 2001, Stamatellos & Whitworth 2004). The density 
structures are well characterized by Bonner-Ebcrt spheres 
(hereafter, BE): pressure bounded, isothermal solutions to 
the equations of hydrostatic equilibrium that do not in- 
clude the effects of magnetic fields or turbulence. A gen- 
eral picture of the initial stages of low-mass star formation 
involves the evolution of a BE sphere from low central den- 
sity to high central density, perhaps retarded by magnetic 
pressure, until a central hydrostatic core develops. If we 
wish to understand the beginning of this process, we must 
identify and study low density PPCs. 

In this paper we present deep submillimeter maps of 
L1498 (§2.1) with radiative transfer modeling which in- 
cludes a more sophisticated treatment of the heating of 
the PPC (§3). The results of dust continuum modeling 
are discussed in §4. A comparison between the derived 



properties from dust continuum observations of PPCs that 
were observed with ISO and SCUBA is discussed in §4.1.2. 
We also present the highest resolution single-dish N 2 H + 
and C3S maps of L1498 (§2.2, §2.3) with analysis of varia- 
tions in the column density, abundance, and velocity (§4.2, 
§4.3). 

2. OBSERVATIONS 

2.1. SCUBA Continuum Observations 
2.1.1. Reduction 

L1498 was observed simultaneously at 850 and 450 /<m 
with SCUBA on the 15 m James Clerk Maxwell Telescope 
during the nights of August 29 and 30, 1998. A total of 
fifty jiggle maps were made towards the source, each map 
consisting of four 64-point jiggle maps with 256s of ON- 
source integration time. The total ON-source integration 
time was 3.55 hours. Previous spectral line maps showed 
that L1498 is extended in a SE to NW direction; therefore, 
the chop angle was set to a constant position angle of 20° 
with a 120" chop throw to chop perpendicular to the ma- 
jor axis of the core. Non-azimuthal chopping results in a 
slight reduction in signal-to-noise in the map since a sym- 
metrical chopping pattern does not lie at constant airmass. 
A preliminary map made in April, 1998 showed that the 
source was detected to the edge of a single 64-point jiggle 
map ficld-of-view (2.3'). Therefore, the map was extended 
using three offset 5-pointing maps, each with 30" spacing 
(see Figure la). The five-pointing maps were centered at 
(0",0"), (+30", -30"), and (-30", +30") from the central 
coordinates of J2000.0 a = A h 10 m 52.5 s , 6 = +25° 9' 
55". The final maps spans ±140" in right ascension and 
declination. 

Each map was reduced using the standard SURF 
(SCUBA User Reduction Facility) reduction routines (Jen- 
ness & Lightfoot 1997). Each 64-point jiggle map was flat- 
fielded and corrected for chop throw, extinction, and sky 
noise. The 450 fim images were reduced using a set inten- 
sity scale to ensure consistent identification of improper 
sky noise subtraction. The telescope pointing was checked 
every hour. The largest pointing shift was 2.5"; therefore, 
the standard pointing offsets from the five-pointing maps 
were used to shift and add together the individual maps. 

Skydips at 850 fim and 450 fim were performed ev- 
ery hour during both nights. These skydips are com- 
pared to the opacity measured every 10 minutes at 225 
GHz (tcso) and 350 /Ltm from tippers located at the 
Caltech Submillimeter Observatory. The tcso data was 
only available for the night of August 30. Using the re- 
lationship derived by Archibald et al. (2002) between 
tcso and the skydip-determined opacity at 850 ^m (rsso), 
t 850 = (3.99 ± 0.02) {tcso - 0.004 ± 0.001), we found ex- 
cellent agreement between our skydips and the scaled 850 
/jm opacity (Figure 2). Therefore, we use the the scaled 
tcso values to correct for the 850 /im opacity in individual 
jiggle maps on the night of August 30. We linearly inter- 
polated between 850 ^m skydips for the night of August 
29. 

Since the sky opacity at 350 /j,m and 450 (im is very 
sensitive to short term variations in atmospheric precip- 
itable water vapor, we used the tipper opacity at 350 /an 
to monitor the variability of T450 between the hourly sky- 
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dips. We found an average ratio of T450 skydips to T350 
of 0.71 ± 0.07 for both nights. The 450 /im opacity was 
corrected for each jiggle map using the scaled T350 opac- 
ity interpolated from the tipper measurements bracketing 
each jiggle map. 

2.1.2. Images 

The reduced maps are shown in Figure 1. L1498 was 
detected at both wavelengths with a peak signal-to-noise 
of 18 and 13 for the 850 and 450 /im maps, respectively. 
The contours are spaced by 3a for both images. The sub- 
millimeter images are characterized by an amazingly flat 
intensity plateau. A 1200 ^im continuum map, smoothed 
to 20" resolution (Tafalla et al. 2004), is also shown for 
comparison. The FWHM intensity contour of the 1200, 
850, and 450 /im maps are very similar and well fitted by 
an ellipse with a major axis (a) of 197", minor axis (6) of 
108", and position angle of 122°. The 850 /im and 450 /im 
ellipses have similar centroids (+11", —16") with an uncer- 
tainty of ±4" with respect to the pointing center (§2.1.1). 
The HWHM radius of the core, defined by the geomet- 
ric mean of the major and minor axes (R1/2 = Vab/2), 
is 73". While this radius may be used an a representa- 
tive size of the core, it is important to realize that the 
FWHM contour is asymmetric with an aspect ratio (a/b) 
of 1.8. Furthermore, the 1200 /im and 850 /im maps are 
non-axisymmetric with a sharper gradient in the intensity 
along the northeast edge compared to the southwest edge. 

The ISO map at 170 /im is also shown in Figure lb 
(Ward- Thompson et al. 2002). This is the shortest wave- 
length at which the far-infrared emission has been ob- 
served to peak on the core; ISO 90 /im observations detect 
a more diffuse emission peak approximately 3' south of the 
submillimeter centroid. The scale of the 170 /im map is 
±300", twice the scale in the (sub)-millimcter maps. The 
ISO beam was very large (« 80") at 170 /im; therefore, 
the map does not show the non-axisymmetric structure 
observed on smaller scales at 850 and 1200 /im. The core 
is elongated, but with a slightly different position angle 
on larger scales (PA ~ 150°). Far- infrared emission is de- 
tected on very large scales, up to 300" (42,000 AU) from 
the centroid of the submillimeter emission. 

Azimuthally-averaged, normalized radial profiles of the 
850, 450, and 1200 /im maps were calculated. The images 
were re-binned to pixels with half beam spacing (7" at 850 
microns and 3.5" at 450 /im) corresponding to the Nyquist 
sampling limit in the map. The centroid was chosen to be 
the average of the 850 and 450 /im centroids (+11", —16"). 
The 1200 /im profile was binned at 10" since the image was 
smoothed by a 20" Gaussian (Tafalla et al. 2004). The flat 
intensity plateau has the same extent in all of the profiles. 

Since L1498 is a very elongated core, azimuthally av- 
eraged profiles will be weighted towards a flatter inten- 
sity profile (see J0rgensen et al. 2002). This effect may 
be quantified by considering the normalized profile within 
sectors centered on the major and minor axes. If the sec- 
tor centered on the major axis is chosen with the opening 
angle given by, 

# maj = 2tan" 1 (&/a) , (2) 
where b/a is the ratio of the minor to major axes (f? m ; n + 
#maj = tt), then the areas of sectors are equal. The 



azimuthally averaged and sector-averaged profile will be 
modeled separately in §3. 

2.1.3. Calibration and the SED 

Uranus and AFGL618 were observed on August 29 and 
30 as flux calibrators and for beam profiles at 850 and 450 
/im. The total flux of Uranus was 74.0 Jy at 850 /im and 
195.5 Jy at 450 /im. The apparent diameter of Uranus was 
3"7 effectively making Uranus a point source at both wave- 
lengths (see Shirley et al. 2000 for an analysis of Uranus' 
finite size on the beam size). The flux of AFGL618 was as- 
sumed to be 4.56 + 0.17 Jy/beam at 850 /im and 11.2 + 1.4 
Jy/beam at 450 /im (SCUBA Secondary Calibrator Web- 
page 4 ). Since we are interested in the total flux observed 
towards L1498, we calibrated the flux in apertures of 40" 
and 120". Flux densities are calculated for an aperture di- 
ameter, 6, using S„(\,6) = C{\,6)V(\,6)e- TSCCZ , where 
V(A, 9) was the voltage measured at wavelength A in an 
aperture of diameter 9 at airmass secz and C(A, 9) is the 
calibration factor. The fluxes for L1498 are included in 
the source properties in Table 1. The uncertainty on the 
calibration factor was calculated by propagating the errors 
using, 



~ 2 C 2 



c 2 



+ 



Run Avg 



/ 2 2 

(sec z <7 r 



Source 



+ 



V 4 2 



C 2y2 

(3) 

where S„ is the flux density (Jy), C is the calibration factor 
(Jy/V), V is the total voltage observed within the aper- 
ture (V), and Ai is the area of a pixel included within the 
aperture. 

The SCUBA beam profiles were determined using ob- 
servations of Uranus at 850 and 450 /im. Uranus was ob- 
servable only at the beginning of our observing shift and 
only at high airmass (sec z <~ 2). The Uranus profile can- 
not be traced beyond 60" at 850 /im and beyond 50" at 
450 /im. Unfortunately, AFGL618 is not strong enough to 
be used for extended beam profiles (6 > 25"). The Uranus 
profiles are very similar to the average Uranus beam pro- 
file determined from maps made in April, 1998 (Shirley 
et al. 2000). Since all of our observation were made dur- 
ing the morning observing shift (1:30AM HST - 9:30AM 
HST), the dish should have cooled sufficiently such that 
the beam shape is stable. No significant trends were seen 
in the April 1998 beam maps observed during the morn- 
ing observing shift; therefore, we feel justified in using the 
Uranus beam profile for all second shift observations. 

The SED was constructed from our SCUBA observa- 
tions and a search of the literature (Table 1). The SED 
includes previous (sub)-millimeter observations by Ward- 
Thompson et al. (1994) as well as ISO far-infrared obser- 
vations at 170 and 200 /im (Ward-Thompson et al. 2002). 
The inclusion of the far- infrared fluxes is extremely impor- 
tant for bracketing the peak of the SED. 

2.2. N 2 H+ Observations 

2.2.1. BEARS 

Lynds 1498 was mapped with BEARS, a 25-beam 
double-sideband SIS receiver array that operates at 3mm 
(Sunada et al. 2000), on April 26, 2004 at the Nobeyama 
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Radio Observatory 45 m telescope. The average beam- 
size within the array is 17'.'8 at 93.17 GHz (Tatematsu et 
al. 2004); however, each beam in the array is separated by 
slightly more than 2# fwhm (4171 ). We observed the source 
using a dithered 4x4 position-switched grid on the sky 
spaced by 20'.'55. The array was rotated by ±90° after each 
dithered map to change relative beam positions with re- 
spect to the source. A digital autocorrelator backend with 
37.8 kHz resolution (w 0.12 km/s) was used. BEARS was 
tuned to the frequency of the strongest hyperfine compo- 
nent (93.1737767 GHz JF X F = 123 -> 012, see Appendix 
A). The resulting N2H + map has 100 spectra separated 
at nearly full beam resolution and represents the highest 
resolution N 2 H+ spectrum obtained towards L1498. 

BEARS tuning was consistent across the array with an 
average T sys of 220 K. The intensity was calibrated for the 
receiver sideband ratio by observing the bright high-mass 
star-forming core, G40.50+2.54, with every beam in the 
BEARS array and the Nobeyama facility single-sideband 
S100 receiver. Spectra from the S100 receiver were cali- 
brated on the TjJ scale using the standard chopper-wheel 
method (Penzias & Burrus 1973). The calibration source 
was observed twice with the array rotated by 90° to min- 
imize systematic calibration errors. The resulting calibra- 
tion is good to approximately 20% . The spectra were 
converted to the scale by assuming the BEARS side- 
lobes are negligible and using i] m k = 0.51 (Tatematsu et 
al. 2004). 

Pointing observations were performed every hour using 
the SiO i> = 1, J = 1 — » mascr transition towards NML 
Tau. The largest pointing shifts were 6" with average cor- 
rections of 4" and 3" in azimuth and elevation respec- 
tively. The weather was clear with winds typically below 
6 m/s during L1498 observations. Conditions were good 
for N 2 H + observations and the absolute pointing errors 
are small compared to the BEARS bcamsize. 

The data were baselined, summed, calibrated, and in- 
terpolated onto a regular spatial grid using the NewStar 5 
software package of the Nobeyama Radio Observatory. 
The reduced maps were written out to the group FITS 
format (Greisen & Harten 1981) and read into AIPS++. 
Fitting of the hyperfine components were performed with 
Glish scripts in AIPS++. The group FITS data cube was 
also written to a series of ASCII files for calculation of 
baseline rms and the integrated intensity. 

2.2.2. CSO Observations 

The N2H + J = 3 — > 2 transition was observed in 
position-switched mode towards the centroid continuum 
position of Lynds 1498 with the Caltech Submillimeter 
Observatory on the night of July 24, 2004. The 230 GHz 
double sideband SIS receiver (Kooi 1992) with a 50 MHz 
AOS backend were used. The average T sys was approxi- 
mately 430 K. The spectral resolution was 147 kHz (0.15 
km/s) after smoothing to the actual resolution of the spec- 
trometer (3 channels). The receiver was tuned to the 
JPL line catalog frequency, 279.511701 GHz, for N 2 H+ 
J = 3 — > 2; this frequency is within 160 kHz of a blend con- 
taining the strongest hyperfine component (279.5118631 
GHz, JFiF = 345^23 4; Luca Dore, priv. comm.). 

5 http://www.nro.ac.jp/ nro45mrt/NEW45m/NewStar/index.html 

6 Frequencies for transitions of C3S and its isotopologues are from 
http://www.phl.uni-koeln.de/vorhersagcn/ 



Pointing and main beam calibration were performed on 
Venus. The main beam efficiency was determined to be 
Vmb = 0.60 ±0.06 at 279 GHz for a 2 6'.' 8 beam. The point- 
ing was consistent with a largest pointing shift of less than 
3" in both azimuth and zenith angle during the observa- 
tions. The spectra were smoothed, baselined, and summed 
using the CLASS software package (Buisson et al. 2002). 
L1498 was observed for a total of 10 minutes ON-source. 
The final spectrum has a baseline rms of cr(T^) = 45 mK. 

2.3. C3S Observations 

The C 3 S J = 4 -> 3 transition at 23.1229820 GHz 6 
was mapped using the 100m Robert C. Byrd Green Bank 
Telescope on the nights of March 14 and 15, 2005. The 
spectrometer was setup to observe 4 IFs with 2 polariza- 
tions at 6.125 kHz resolution (0.079 km/s) using a single 
beam of the K-band receiver. The four IFs were centered 
on the C 3 S C 3 34 S, 13 CCCS, and C 13 CCS J = 4 -» 3 
transitions. The observations were frequency switched by 
+4 MHz at a frequency of 4 Hz. Initially, a pointed map 
was made using a 7 x 5 grid, with 60s integration time 
per pointing, spaced at 33" (full beam spacing) centered 
on the SCUBA dust continuum centroid position. Three 
smaller 3x3 maps with 16'.' 5 spacing were centered at 
the offsets (0,0), (+49'.'5,-16'.'5), and (-4975, +1675) with 
respect to the dust continuum centroid. 

The raw frequency-switched spectra were folded and cal- 
ibrated on the Ta scale using the DISH software package 
within AIPS++. The data were then written to ASCII 
files where opacity corrections were applied (T A ) and RR 
and LL polarizations were summed with 1 / cr^ aseline weight- 
ing. The opacity at 23.1 GHz was determined using 
the weather model of R. Maddalena (private communica- 
tion) which averages the precipitablc H 2 measurements 
from three nearby towns (Elkins, Hot Springs, and Lewis- 
burg, West Virginia). The average zenith opacity was 
T23.1 = 0.049±0.002 on March 14 and t 22 . 5 = 0.106±0.002 
on March 15 during the observations. Final conversion to 
the T m b scale was achieved using position-switched scans 
of the quasars 3C48 and 3C286. Unfortunately, the aper- 
ture efficiency was not properly determined at 23.1 GHz 
on March 14; therefore, we have linearly interpolated from 
measurements made at 9 GHz (Langston et al. 2004) 
and 46 GHz (March 18, 2005) to find r/ a = 0.54 ± 0.11. 
The GBT has an active surface which mitigates varia- 
tions in the aperture efficiency with elevation and tem- 
perature. We have applied a 20% systematic uncertainty 
to account for these effects. The main beam efficiency 
is then r] m b = 0.74 ± 0.15. The T m b temperature scale 
may be compared with observations on the scale if the 
sidelobe pickup is negligible (see Chapter 8.2.5 of Rohlfs 
& Wilson 2000); this is a reasonable assumption for the 
K-band receiver on the GBT. 

Pointing and focusing were performed once per hour. 
The sky was clear and the wind was less than 4 m/s on 
both March 14 and March 15 during the observations. 
Since the observations were performed at nighttime, the 
temperature variation was small (AT < 1.9°C). The point- 
ing and focus were very stable with a largest pointing shift 

Cologne Database for Molecular Spectroscopy (Miiller et al. 2001): 
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of 7". 

3. MODELS 

3.1. ID Continuum Radiative Transfer 

The dust continuum observations arc modeled using 
the one-dimensional radiative transfer code CSDUST3 
(Egan, Leung, & Spagna 1988) that self-consistently calcu- 
lates the temperature distribution for a given density dis- 
tribution, dust opacity, and ISRF. CSDUST3 simulates 
anisotropic scattering which may affect the temperature 
distribution in the outer regions of the core, a feature 
that is not present in the radiative transfer code of Ivezic, 
Nenkova, & Elitzur (1999). Once the temperature distri- 
bution has been calculated, the intensity profiles at 850 
and 450 /im are created by convolving the model inten- 
sity profiles with the measured beamshape and simulating 
chopping. The model SED is also calculated and com- 
pared to the observed fluxes. Simultaneous modeling of 
the intensity profiles and SED result in nearly orthogo- 
nal constraints on the shape of the density profile and the 
total mass for a given opacity choice (see Shirley et al. 
2002). The modeling procedure for PPCs is described in 
more detail in Evans et al. (2001). 

We shall use a family of BE spheres (Ebert 1955, Bonnor 
1956) to attempt to fit the L1498 observations since BE 
spheres were used by Evans et al. (2001) to successfully 
fit submm continuum observations of other PPCs. BE 
sphere are solutions to the pressure bounded equations of 
hydrostatic equilibrium that ignore the effects of magnetic 
fields and turbulence. The BE density profile is flat in 
the inner regions and asymptotically approaches r~ 2 in 
the outer regions. The size of the flat density plateau de- 
creases as the central density of the BE sphere increases 
(density FWHM = R(\n c ) ~ 1/y/fiZ, Chandrasekar & 
Wares 1949). The temperature structures are character- 
ized by decreasing dust temperature from outside to in- 
side due to extinction of heating radiation from the ISRF. 
The temperature gradient from outside to inside increases 
for more centrally condensed BE spheres. Isothermal BE 
spheres may be parametrized by a single variable, the cen- 
tral density (n c ) for a fixed temperature (10 K). Evans et 
al. (2001) showed that non-isothermal correction to the 
shape of the BE profile was negligible for profiles used to 
describe PPCs. Isothermal BE spheres are calculated with 
central densities ranging from logn c = 3.5 to logn c = 7.0. 
Since L1498 is embedded within a larger molecular cloud 
with Ay « 1-3 mag (Cambresy 1999), the power- law por- 
tion of the BE profile is not continued below 10 3 cm~ 3 . 

Since PPCs lack an internal source, the heating is from 
the ISRF. The spectrum of the ISRF is derived from 
COBE data (Black 1994) with ultraviolet wavelengths 
constructed using an empirical description of the radia- 
tion field that reproduces the ISRF of Draine (1978; see 
van Dishoeck 1988). This spectrum is the same as used in 
the models of Evans ct al. (2001). The shape of the ISRF 
spectrum may be modified by two effects: the intensity of 
the ISRF (UV to far-infrared) varies at different locations 
throughout the Galaxy; and the short wavelength inten- 
sity of the ISRF (UV to near-infrared) may be modified by 
dust extinction (see Figure 3) . The scaling of the intensity 
is parametrized by s; sr f while the amount of dust extinction 
at the core's outer radius is parametrized by Ay. We use 



dust opacities appropriate for the general ISM (graphite 
and silicate based grains, Draine & Lee 1984) to modify 
the ISRF spectrum. The combination of scaling the ISRF 
and extinction at the outer core radius is a more realistic 
description of the incident intensity than has been used in 
previous models of PPCs. 

The reprocessing of the ISRF is strongly affected by the 
dust opacity within the core. We model the emission us- 
ing eight different opacities models (see Table 2 and Figure 
3). The opacities of Ossenkopf & Henning (1994) are de- 
rived from models of grains that have coagulated for 10 5 
years and that have acquired varied depths of ice via gas 
adsorption: OH2 no ice mantles (column 2 Table 5 of Os- 
senkopf & Henning 1994), OH5 thin ice mantles (column 5, 
Table 5), and OH8 thick ice mantles (column 8, Table 5). 
The opacities of Mathis, Mezger, & Panagia (1983) are de- 
rived from an empirical fit to dust properties. We also test 
the opacities of Weingartner & Draine (2001) which de- 
scribe silicate and graphite grains, including a population 
of small grains and PAHs, appropriate for the ISM. The 
opacities of Pollack (1994) describe an alternative grain 
composition based on silicate grains (olivine and orthopy- 
roxene), iron compounds (troilite and metallic iron), and 
various organic C compounds. The mass opacity at 850 
/im, k„(850), varies by more than one order of magnitude 
between the opacity models (Table 2). A constant opacity 
with radius (R e [0, R a ], R Q is the outer core radius) is 
assumed throughout the core. 

3.2. L1498 Model 

A total of 770 radiative transfer models were run that 
vary n c , R 0l s; sr f , Ay, and n v . For each model, the reduced 
X% given by 

a _f (^-af)'/^) 2 

Xr TV ' [ ' 

i=l 

was calculated for the 850 and 450 normalized intensity 
profiles and for the SED. A \ 2 is also calculated for the flux 
in a 120" aperture at 850 /im. The best-fitted models are 
determined by finding the intersection of the sets of Xr f° r 
each quantity within the typical range Xr € [0,1]. Since 
the errorbars of the intensity profiles include azimuthal 
variations as well as intensity uncertainties, and since the 
errorbars of the SED include an estimate of systematic 
uncertainty, it is possible for Xr < 1- 

The first grid of models varies n c and R Q for s; sr f = 1, 
Ay = 1 , and OH5 opacities (see top row of Figure 4) . The 
best-fitted model to the azimuthally averaged, normalized 
intensity profiles is n c = 1 x 10 4 cm~ 3 . If we re-calculate 
the Xr(^85o) using the sector-averaged, normalized inten- 
sity profiles, then a wider range of central densities fit, 
n c = lt ; 5 x 10 4 cm~ 3 (see dashed profiles in Figure 5d). 
The best fitted central density does not strongly depend 
on the size of the outer radius; radii from 30,000 AU to 
60,000 AU are well fit. Since emission is clearly detected 
100" from the center of the SCUBA maps, R Q must be 
greater than 15, 000 AU. Furthermore, the ISO 170 /im 
image (Ward-Thompson et al. 2002) detects dust emis- 
sion at radii greater than 200" (28, 000 AU) from the cen- 
ter; therefore, the best-fitted model R are consistent with 
submillimeter and far-infrared observations. 

Since the SCUBA continuum centroid and 1.2 mm con- 
tinuum peaks are different by 31", we have also compared 
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the dust model intensity profiles to azimuthally averaged 
intensity profiles centered at the 1.2 mm peak. The best- 
fitted model are n c = 1 — 3 x 10 4 cm~ 3 Bonnor-Ebert 
spheres. This agreement with the best-fitted models to- 
ward the SCUBA centroid position is a result of L1498's 
very flat submillimeter intensity profile that does not vary 
strongly with different central positions within 30" of 
the SCUBA continuum centroid (see Figure 1). 

The second grid of models varies s; sr f, and Ay for the 
n^lx 10 4 BE sphere, R Q = 20, 000, 40, 000, and 60, 000 
AU, and the different opacity models. The Sj sr f and Ay are 
varied on non-linear grids with values between [0.2, 10] and 
[0.1, 10] respectively. The middle and bottom rows of Fig- 
ure 4 display xf(SED) and Xr(^85o) f° r OH5 dust opaci- 
ties. Changing Si sr f and Ay has a negligible effect on the 
best fitted normalized intensity profiles; the n c = 1 x 10 4 
BE sphere is the best fit for all opacity models. Previ- 
ous dust continuum models relied only upon submillime- 
ter fluxes to constrain the heating properties of the ISRF. 
However, Xr(S&5o) alone do no strongly constrain Sj sr f and 
Ay (bottom row, Figure 5); a wide range of values provide 
good fits. When we include the complete SED, the s; sr f is 
constrained to be < 2 for the best-fitted models (middle 
row, Figure 5). 

The only dust opacity models that provide satisfactory 
fits to 7s50, $850, and the SED are the Ossenkopf & Pfen- 
ning opacities. None of the model xl{SED) for WD3, 
WD4, Pollackl, or Drainc-Lcc opacities are below 6. These 
models underestimate the fluxes observed at far-infrared 
through millimeter wavelengths for the 1 x 10 4 cm~ 3 BE 
sphere. This is not surprising since the mass opacity at 850 
/im for these models is 3 to 8 times smaller than the mass 
opacity of OH8 dust. Furthermore, Ossenkopf & Hcnning 
opacities for grains with ice mantles (OH5 and OH8) pro- 
vide better fits to Sgso and the SED than grains without 
ice mantles (OH2). The OH2 opacities tend to overesti- 
mate 5*850 while underestimating the far-infrared fluxes, ft 
is encouraging that the dust opacities for grains with ice 
mantles provide the best fits since these models are qual- 
itatively consistent with a physical model of L1498 that 
accounts for the observed gas depletion of species such as 
CO, H 2 CO, CS, and CCS. 

The best-fitted models from all of the model grids are 
listed in Table 3 and shown in Figure 5. The models fall 
into two categories: OH5 opacities with Si sr f G [1,2] and 
Ay g [4, 8]; and OH8 and OH5 opacities with s isr f g [0.5, 1] 
and Ay g [0.5, 2]. The fluxes and intensity profiles are well 
fit at all wavelengths except 1.2 mm for all of the models 
All of the best-fitted models have an outer radius of 60, 000 
AU (0.29 pc). Since a large outer radius is preferred, it 
seems unlikely that the high Ay models are physical. An 
optical extinction map of the Taurus molecular cloud lim- 
its the large-scale extinction around L1498 to Ay < 3 mag 
(Cambresy 1999). Indeed, for Ay = 6, the column density 
of material outside R a would be Nh 2 ~ 1 x 10 22 cm~ 2 ; 
this value is larger than the column density derived from 
a 1 x 10 4 cm" 3 BE sphere with R a = 60, 000 AU. We prefer 
the low Ay models as the best, most physical fit. 

All of the best-fitted models are characterized by dust 
temperatures between 10 and 11.1 K (Figure 5a). The low 
Ay OH5 and OH8 models have temperature gradients of 
0.3 < AT[ Q Ro ] < 0.5 K. These are nearly isothermal tem- 



perature profiles compared to previous models of PPCs 
(AT [0;i jj > 1 K, e.g. Evans et al. 2001). The low density 
of L1498 combined with the inclusion of dust attenuation 
at the outer radius results in the nearly isothermal tem- 
perature profiles. 

4. ANALYSIS 

4.1. Derived Properties from Continuum Observations 
4.1.1. L1498 

The bolometric luminosity was calculated by integrating 
under the spectral energy distribution of L1498 (Table 1) 
and is 0.12 ± 0.02 L©, low for a PPC, but not the lowest 
reported luminosity (see Table 4). The bolometric tem- 
perature, the temperature of a blackbody with the same 
mean frequency as the observed SED (Chen et al. 1995), 
is 13.9 ±3.1 K. This is one of the lowest bolometric tem- 
peratures reported for a PPC (see Table 4). 

We calculate the mass from the flux at 850 /im in a 120" 
aperture and assuming a constant density, isothermal dust 
temperature, and optically thin emission using, 



exp iw) - 1 



(5) 



Based on ID dust radiative transfer models, a dust tem- 
perature of 10.5 ± 0.5 K is used for L1498. The resulting 
dust mass, Mn(10.5K) = 0.67 ± 0.16 M© for an aver- 
age of OH5 and OH8 opacities (7c7;(850) = 0.020 ± 0.002 
cm 2 g -1 ). The errors in Md are propagated from S u , T d , 
and 7^7(850). L1498's dust mass is consistent with and 
slightly lower than previously observed pre-protostellar 
cores ((M D (f 0.5)) = 0.9 ± 0.3 M , Table 4). 

We also calculate the aperture-averaged column density 
from the 850 /im fluxes using 



Nh 2 



2S> 2 



exp 



/ hv 
\kT~d 



1 



(6) 



hv 3 fimnKvivOap 
The column density at the SCUBA continuum centroid 
is (1.25 ± 0.71) x 10 22 cm~ 3 in a 17"6 aperture, where 
errors are propagated in S v , T4, and k7J\ This aperture 
size matches the BEARS beamsize allowing a direct com- 
parison between dust and N 2 H+ column density (§4.1.2). 
The column density corresponds to a visual extinction of 
Ay = 13 mag (Whittet 2003). This extinction is compa- 
rable to the peak extinction observed using near-infrared 
extinction mapping towards the low-density PPC, Coal- 
sack G2 (Ay = If .5 mag, Lada et al. 2004). Modeling of 
the Coalsack G2 extinction profile indicates a peak density 
of w 10 4 cm -3 , corroborating the low density derived for 
L1498 from our radiative transfer models. 

Dust opacities tend to follow a single power-law with 
frequency at submillimeter wavelengths (k v cx v@, A g 
[350,1300] /im). The observed dust emissivity index, (3, 
for a core observed at 850 and 450 /im is given by 



log 



5*450 



log 



' 850 > 
v450/ 



3+- 



log 



cxp(hug 5a /kT d )-l 

CXp(ft^ 4 5o/feT ti )-l 



log( 



850 1 
15U J 



(7) 

where the term in square brackets is the spectral index 
between 450 and 850 /im , 0450/850 1 and the term in curly 
braces is the correction for failure of the Rayleigh- Jeans 
approximation. Assuming a dust temperature of 10.5 ±0.5 
K, we find a£§g /850 = 2.91 ±0.61 and (3 = 2.44±0.62. The 
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spectral index is higher than all of the PPCs observed by 
Shirley ct al. (2000) except for the two cores associated 
with L1689A. The (3 is larger than the theoretical (3 for 
dust opacities used in modeling but is consistent, within 
errors, with all of the dust opacity models except OH2 
and WD4 (see Table 2). While OH5 and OH8 opacities 
give good fits, within errorbars, to the SED, they consis- 
tently underestimate the flux at 450 /xm while matching 
the 850 /im flux. Submillimeter measurements of a and 
P do not discriminate well between opacity models; ra- 
diative transfer modeling of the complete SED is a better 
discriminator. 

4.1.2. Comparison with Nearby PPCs 

The properties of Lf498 may be placed in context by 
comparing to other PPCs that have been observed with 
ISO and SCUBA. Table 4 lists observed and derived prop- 
erties of 10 nearby (D < 200 pc) sources. Including the 
170 /im ISOPHOT flux gives at least one flux point on 
the Wien side of the SED, resulting in a more accurate 
calculation of L bo i- Each source SED contains at least 3 
flux points (170, 200, and 850 /im). We also analyzed the 
SCUBA 850 ^m images for each source. All of the SCUBA 
observations are published except for three cores, L183, 
L1709A, and L63 which were obtained from the CADC 
SCUBA archive. The reduction procedure for these cores 
is identical to that described in §2. We have calculated 
the size, aspect ratio, bolometric luminosity, standard 
evolutionary indicators (T bo i and r Lb ° l ), and the dust- 
determined mass from the 850 /im flux. This is a small 
sample but probably representative of nearby PPCs. The 
CADC SCUBA archive contains many more observations 
of PPCs; however, far-infrared observations are needed to 
calculate L bo i- 

L1498 is the largest nearby PPC in the sample by 50% 
and nearly twice as large as the average core size. L1498 
also has one of the lowest masses (0.67 M ) of the sample. 
This is consistent with the interpretation that L1498 is a 
low density core with a uniquely flat intensity profile. In 
contrast, the average nearby PPC is a core with 0.14 L Q , 
0.8 M©, and elongated ((a/b) = 1.6) with a mean radius 
of 6000 AU. 

PPCs are clearly separated from Class sources on a 
Tbci vs. 773- diagram (Figure 6). All of the PPCs have 

< 25 whereas the observed Class 



L ,s n 



low T bol < 20 K and 
sources all have T boi > 20 K. T bo i and for PPCs 

do not correlate. This is not surprising since these evolu- 
tionary indicators are strongly influenced by the strength 
of the ISRF; this is an environmental effect not directly 
related to the evolutionary state of the core. A better 
scheme for characterizing the evolutionary state of PPCs 
would be the comparison of the dynamical maturity deter- 
mined from radiative transfer models (n c ) and the chemi- 
cal state of molecules that are abundant at different times. 
For instance, comparisons of [CCS]/[N 2 H + ], comparisons 
of [N 2 D+]/[N 2 H+], and the depletion fraction of, fo, of 
C^O and H 2 CO are promising candidates. The evolu- 
tionary state of L1498 is discussed in §4.3. 

4.2. Derived Properties from N%H + Observations 

The seven hyperfine lines of N 2 H + J = 1 — > were well 
detected over an extent of ±80" in the map (Figure If). 



The peak T* A in the map is coincident with the contin- 
uum ccntroid (Figure 7a). We can estimate the effective 
density needed to excite an easily observable 1 K N 2 H + 
line (see Table 1 of Evans 1999). The critical density of 
the J = 1 — > is n cr u = 3.8 x 10 5 cm~ 3 (collision rates 
determined from Turner 1995). By analogy with HCO + , 
the resulting effective density (corrected for optically thin 
hyperfine structure) is n e // « 2x 10 4 cm -3 , consistent 
with the derived central density. 

The integrated intensity is calculated by integrating 
over the seven hyperfine components. The I(T^) at the 
continuum centroid is also the peak in the N 2 H + map, 
J(T*) = 3.92±0.48 K km/s. The FWHM integrated inten- 
sity contour agrees very well with the FWHM continuum 
contour in size and position angle. The N 2 H + emission 
appears to trace the general shape of the continuum emis- 
sion, but does not peak as sharply along the north-eastern 
ridge as is seen in the 850 ^m and, especially, in the 1200 
/zm images. 

The hyperfine structure was fitted with seven Gaussian 
functions, constrained by the hyperfine transition veloci- 
ties (Appendix A), using the profile fitter task in AIPS++. 
The average linewidth was calculated to be 0.22 ± 0.02 
km/s over the map using the isolated hyperfine compo- 
nent {JFxF = 101 -> 012). This linewidth is slightly 
smaller than the linewidth of 0.28 km/s found by Caselli 
et al. (2002) using the FCRAO. The non-thermal con- 
tribution to the linewidth is Av nt — 0.13 km/s. A plot 
of the linewidth versus radius (Figure 7b) indicates that 
the linewidth does not vary significantly out to 15, 000 AU 
(w 110"). Therefore, the average gas motions in L1498 
have low turbulence and are similar throughout the core. 
This result agrees with the Av-r relationship towards other 
PPCs in Taurus observed with BEARS (Tatematsu et al. 
2004). 

We calculate the virial mass using the Doppler linewidth 
and assuming an isothermal dust temperature, 



(6R 



FLAT 



9R 1/2 ) 



8Gln2 



+ 



Av 2 



fimn 8 In 2 



HI a 



(8) 

For a more accurate calculation, we split the radial profile 
into two sections. The first corresponds to R < Rflat 
(the radius where n — 0.9n c ) which is approximated by a 
constant density. For a 1 x 10 4 cm~ 3 BE sphere, Rflat is 
520 AU. Within the second section, Rflat < R < R1/2, 
the density can be approximated by an inverse power law 
with r~ 2 (R 1/2 = 10200 AU, §2.1.2). The term in square 
brackets converts the observed N 2 H + linewidth to a 3D 
velocity dispersion (Shirley et al. 2002). Using the Av 
determined above, we find M„j r = 0.84 ± 0.21 M Q with 
errors propagated in i?i/ 2 , Ti SO , and Aw. This mass agrees 
with the mass determined from dust continuum emission 
(§4.1.1) indicating that L1498 is likely a bound core and 
that the best-fitted Ossenkopf & Hcnning dust opacities 
are appropriate. 

The N 2 H+ J = 1 — > transition is observed to 
have moderate optical depths towards low-mass cores (e.g. 
Tatematsu et al. 2004, Crapsi et al. 2005); therefore, 
we must calculate the optical depth at each position in 
the map to derive the column density. For a position- 
switched observations, theoretically we observe the quan- 
tity T R = {J v {T ex )~ J„(T 6ff ))[l-exp(-T)], where J V {T) = 
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{hv/k)/[exp(hv/kT) — 1]. This Tr is related to the actual 
observed quantity, through source-beam coupling and 
through the efficiency of the antenna (including spillover 
and atmospheric attenuation, see Kutner & Ulich 1981). 
Since the hyperfinc lines lie close together in frequency 
(Av < 4.65 MHz), we assume that the proportionalities 
for each line are constant (Keto et al. 2004). Therefore, 
we find the optical depth in the strongest hyperfine com- 
ponent (r m , JFxF = 1 2 3 — » 1 1 2) by taking ratios of the 
observed antenna temperatures of the i th hyperfine com- 
ponent to the strongest hyperfine component 

where Ri is the ratio of relative strengths of the hyperfine 
lines to the strongest line. Ri is calculated in Appendix 
A. 

The average optical depth in the strongest hyperfine 
line, (r m ), was calculated using four of the six observed 
line ratios and Equation 9. Two of the lines (JF\F = 
112 — > 012 and 110 — ► Oil) consistently gave lower 
(r m ) and higher (r m ) by factors of 0.4 and 2.0 respec- 
tively. This may be due to anomalous hyperfine excita- 
tion, an effect that has been observed towards other PPCs 
(Caselli et al. 1995, Turner 2001). An optical depth was 
determined for 47 positions with suitable signal-to-noise. 
The optical depth is not affected by line blending since 
the linewidth is 0.22 km/s (see Appendix A). Since, each 
hyperfine line has a different optical depth, then correct- 
ing the column density using only r m would overestimate 
the column density. We define the average tau by averag- 
ing the theoretical hyperfine ratios for the seven observed 
transitions (t) = |^JJ,T m = %r m . The average opti- 



(X N2 H+) map = (4.0 ± 3.7) x 10~ 10 relative to H 2 . This 
column density agrees very well with the early-time abun- 
dances of the coupled dynamical-chemical models of Lee 
et al. (2004) for quasi-statically evolving PPCs (§4.3). 
The abundance show no significant variations with radius 
out to 15,000 AU, although the errors on the abundance 
are sizable (Figure 7d). Radiative transfer models of the 
N 2 H + emission should assume a constant abundance of 
w 4 x 10- 10 . 

The N 2 H + J = 3 — > 2 transition was not detected 
towards the continuum centroid position with a lcr(T^) 
baseline rms of 45 mK (Figure 7c); however, Crapsi et al. 
(2005) report a detection with the IRAM 30 m toward the 
1.2 mm continuum peak position (TJ = 380 ± 53 mK; see 
inset in Figure 7c). This is puzzling since the J = 1 — > 
BEARS map shows no enhancement toward that position. 
It is possible that the 1.2 mm continuum peak represents 
the true density peak of L1498. A high-resolution, over- 
sampled map of the N 2 H + J = 3 — > 2 emission is needed to 
resolve this discrepancy. The exact location of the density 
peak will be important for molecular line radiative trans- 
fer modeling, but does not affect the conclusions from dust 
continuum modeling (§3.2). 

The LI 498 N 2 H + observations will be modeled, along 
with other PPCs, in a future paper. 



cal depth in the map is (r) 



map 



0.83 ± 0.54, indicating 



that the N 2 H+ J = 1 — > transition is moderately op- 
tically thick over most of the map. As an example, the 
SCUBA continuum centroid position has an optical depth 
of (r) = 1.11 ± 0.61, resulting in a correction for optical 



depth to the column density of (t)/(1- 



1.66 



+0.43 

-0.38- 



The total column density for N 2 H + is (sec Goldsmith & 
Langer 1999) 7 



N 2 H+ 
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4.3. Derived Properties from C3S Observations 

The C3S J = 4 — > 3 line was strongly detected toward 
L1498 (Figure 8) and the emission follows the same double- 
peaked pattern as previous observations of CS and CCS 
(Figure 1 of Willacy, Langer, & Vclusamy 1998). The 
line intensity is strongest in the southeastern peak (T^ 
= 1.64 ± 0.35 K) at J2000.0 a = 4 h 10 m 56.9 s , 5 = +25° 
9' 22'.'5. Since the C3S emission has an intensity saddle at 
the SCUBA dust continuum centroid, it is another excel- 
lent example of strong chemical differentiation in L1498. 
The southeastern C3S peak does not exactly correspond 
to the Willacy, Langer, & Velusamy CS or CCS peak, but 
instead lies between them and closer to the CS peak. Com- 
bination of the GBT observations with higher resolution 
observations (e.g., VLA) would improve the peak position 
and may confirm slight chemical differentiation between 
all of the observed sulfur-carbon chain molecules. 

The C3S J = 4 — > 3 linewidth was measured using the 
Gaussian fitting routine is AIPS++ toward 42 positions 
in the map with detections greater than 3a. The average 
linewidth is Av = 0.21 ± 0.03 km/s, consistent with the 
N 2 H + linewidth. No significant variation of Av is also seen 
for C3S with radius (Figure 8b). 

The column density of C3S was estimated using the 
same method as for N 2 H + (Equation 10). We estimate 
the optical depth of the C3S J = 4 — > 3 with 10 minute 
observations at the southeastern peak position of three iso- 
topologues. The detection of C3 34 S J = 4 — > 3 (Figure 8c; 
T^(C 3 S)/T^(C3 34 S) = 19.7±3.1) and lack of a significant 
detection (> 3a) of 13 CCCS or C 13 CCS J = 4 -» 3, in- 
dicates, for the standard ISM isotope ratios of [ 32 S]/[ 34 S] 
= 22 and [ 12 C]/[ 13 C] = 77 (Wilson & Rood 1994), that the 
C3S emission is consistent with being optically thin. The 

7 This derivation of the correction for an optically thick column density should not be confused with estimates based on the total optical depth, 



rimh 

(10) 

where Q(T ex ) is the rotational partition function and 
/i = 3.4 D (Green et al. 1974) is the dipole moment. We 
assume an excitation temperature of 10 ± 5 K. Errors in 
A r NaH + are propagated for T ex , (r), and I(T£). The peak 
column density is (8.3 ± 3.4) x 10 12 cm -2 with an average 
column density of (3.5 ± 2.1) x 10 12 cm~ 2 . The peak col- 
umn density agrees remarkably well with the calculation 
of Caselli et al. (2002a) for lower resolution FCRAO 15 m 
N 2 H+ observations towards L1498 ((8 ± 4) x 10 12 cm" 2 ) 
and the intermediate resolution IRAM 30 m observations 
of Crapsi et al al. (2005; (7.1 ± 0.7) x 10 12 cm" 2 ). 

We calculate the abundance, A NaH +(i?) = 
A r N; , H +(i?)/A^H 2 (i?), by comparing the N 2 H+ column 
density to the beam averaged column density deter- 
mined from the 850 dust continuum map (Equa- 
tion 6). The average column density in the map is 



Ttot 



-T m (e.g. Caselli et al. 2002b). 
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peak column density is Nc 3 s = (5.0± 1.8) x 10 12 cm 2 and 
is less than a factor of 2 less than the peak CCS column 
density determined by Wolkovitch et al. (1997) (8.7 x 10 12 
cm -2 ). 

The abundance of C3S was calculated by comparing the 
C3S column density to the column density derived from 
the 850 /xm map in 33" apertures. Since the C3S emis- 
sion is very asymmetrical, we plot the abundance versus 
radius in the major-axis sector (§2.1.2) in Figure 8d. Neg- 
ative radii refer to positions SE of the dust continuum 
centroid. The depletion of C3S is clearly visible with en- 
hanced abundances by factors of 4.8 and 2.3 at the two 
C3S intensity peaks. The abundance is largest in the SE 
peak (X CsS = (4.5 ± 2.4) x 10~ 10 ). The abundance peaks 
are symmetrically located at a radius of 9500 AU from the 
dust continuum centroid position. This radius set a limit 
for the depletion radius of C3S, corresponding to a density 
of 0.8- 1.5 x 10 4 cm~ 3 for the best-fitted BE spheres. Since 
the abundance pattern of C3S is similar to smaller sulfur- 
carbon chain molecules (e.g., CS and C2S), the southeast- 
ern C3S peak of LI 498 would be an excellent site to search 
for larger sulfur-carbon chain molecules such as C4S and 
C 5 S (Gordon et al. 2001). 

4.4. Interpretation of Density Structure 

Radiative transfer modeling of submillimetcr images 
combined with far-infrared and millimeter photometry in- 
dicate that L1498 is a low density, nearly isothermal PPC, 
with a low strength of the ISRF and dust opacities ap- 
propriate for coagulated, icy grains. The central density 
of 1 — 3 x 10 4 cm" 3 is the among lowest central densities 
reported for PPCs, similar to Coalsack G2 (Lada et al. 
2004). This density range agrees with the densities deter- 
mined by Langer & Willacy (2001) from analysis of ISO 
observations of L1498 (1.2 - 5.5 x 10 4 cm -3 with lower 
densities more appropriate for larger core radii). How- 
ever, our modeling results differ substantially from the 
analytical model of Tafalla et al. (2004) which fits the 
radial profile of L1498 at 1.2 mm with a function of the 
form n{r) = n c /[l + (r/i? 1 / 2 ) p ] and assumes an isother- 
mal dust temperature profile of Td(r) = 10 K. Tafalla et 
al. find a central density of 9.5 x 10 4 cm~ 3 , three times 
larger than the maximal BE central density in our mod- 
els. There are a few discrepancies that may account for 
this difference. Tafalla assume a dust opacity of 5 x 10 -3 
cm 2 g _1 at 1.2 mm that is a factor of 2 smaller than OH5 
opacities (1.02 x 10~ 3 cm 2 g _1 ). Since the column density 
is inversely proportional to the opacity, the Tafalla et al. 
opacity would result in a factor of 2 larger central den- 
sity. This does not account for the full difference between 
the results. The Tafalla et al. centroid position is located 
at the peak of the non-axisymmetry which is much more 
pronounced at 1.2 mm than at submillimetcr wavelengths. 
Tafalla et al. calculate a column density of 3 - 4 x 10 22 
cm~ 2 , a factor of 3 higher than our column density de- 
termined at 850 /im. For all of the best-fitted radiative 
transfer models, the 1.2 mm flux is always underestimated 
while the other six fluxes (170, 200, 450, 850, 1100, and 
1300 /urn) are fit within the errorbars. The 1.2 mm flux 
point appears to be anomalously high. It is possible that 
calibration uncertainties may account for the remaining 
difference. 



The low density implies that L1498 is dynamically 
young. Indeed, L1498 may be in a stable hydrostatic state. 
One measure of stability of BE spheres is the density gra- 
dient between inside and outside of the core. If the density 
contrast exceeds 14.3, then the BE sphere is in an unsta- 
ble equilibrium (see Foster & Chevalier 1993). A lower 
density limit of 10 3 cm -3 was assumed for the core in the 
radiative transfer models. This lower limit was chosen to 
be larger than the average density in a molecular cloud 
(10 2 cm~ 3 ) since L1498 is situated within an extended 
condensation of 1 to 3 Ay in the Taurus- Auriga molecular 
cloud (see Cambresy 1999). The average solution (1 x 10 4 
cm~ 3 ) is stable, but this results depends sensitively on the 
density at the outer radius. For the maximal central den- 
sity (3 x 10 4 cm -3 ), the density contrast is 30 and the core 
is unstable. Unfortunately, we cannot strongly constrain 
the density contrast since the submillimeter observations 
are insensitive to structure outside the chop distance (120" 
or 16,800 AU). 

A more quantitative way to test the stability is to 
directly compare the Jeans mass with the mass of the 
BE sphere. The Jeans Mass for T = 10.5 K and n = 
10 4 cm- 3 is Mj w 18T 3 / 2 n-^ 2 = 5.7 M (Spitzer 
1978), while the 1 x 10 4 cm~ 3 BE sphere has a mass of 
M BE « 1.15(T/10K)(P o /10 5 Kcm- 3 )- 1 /2 = 3.6 M Q (Bon- 
nor 1956). The BE mass is within a factor of 2 of the Jeans 
mass, but M BE < Mj. If we use the 3 x 10 4 cm" 3 BE 
sphere, then the mass is comparable to the Jeans mass 
Mj = 3.3 M Q . This mass comparison is consistent with 
L1498 being a core that is marginally stable to gravita- 
tional collapse. 

While L1498 may be dynamically unevolved, it appears 
to be chemically evolved. Species such as H2CO, CS, CCS, 
and C 3 S (e.g., Willacy et al. 1998, Tafalla et al. 2004, 
Young et al. 2004, §4.3) are depleted towards the center 
of the core. We can estimate the depletion timescale from 
the rate equation given in Rawlings et al. (1992). For a de- 
pletion fraction, /d, the adsorption timescale of a neutral 
molecular species, i, is given by, 

1.38 x 10 9 ln(l/,f D K 1/2 

tadsorpW = nS.T 1 / 2 yrS '' ( ) 

where m, is the molecular mass in a.m.u. and Si is the 
sticking coefficient (also see Charnley, Rodgcrs, & Ehren- 
freund 2001). For example, H 2 CO has been observed to be 
severely depleted towards L1498 over a large spatial extent 
(Wang 1994; Young et al. 2004). The adsorption timescale 
is *ad S orp(H 2 CO) w 5.5 x 10 5 years for / D = 1/10, # = 1, 
and T = 10.5 K in the absence of other formation or de- 
struction mechanism. In reality, many chemical effects 
may affect this timescale. Two of the first molecules to 
be depleted are CCS and C3S. Comparisons with detailed 
coupled dynamical-chemical models indicate that the CCS 
becomes depleted on timescales of a few 10 5 years (e.g., Li 
ct al. 2002, Aikawa et al. 2003, Lee et al. 2004). 

A chemical age of a few 10 5 years is consistent with the 
observed abundance of N 2 H + (Aikawa et al. 2003, Lee et 
al. 2004). N 2 H+ is considered a "late time" species since 
it is primarily destroyed in the gas phase by CO and its 
abundance increases after CO is significantly adsorbed (al- 
though dissociative recombination may also be important 
at low densities, see Geppert et al. 2004). Detailed chem- 
ical models indicate that this abundance increase occurs 
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after w 10 5 years (e.g. Aikawa et al. 2003). Clearly, L1498 
must be have been in a stable or very slowly evolving con- 
figuration for at least 10 5 years to display the extreme 
molecular depletion and modest N 2 H + abundance that is 
observed. 

The dust opacities that best-fitted the L1498 SED were 
Ossenkopf & Hcnning opacities (1994) for grains that had 
coagulated for 10 5 years and accreted ice mantles (OH5 
& OH8). These opacities have also successfully fit the 
SEDs of low-mass embedded protostars (e.g., Shirley et al. 
2002) and high-mass embedded protostars (e.g., Mueller et 
al. 2002). The OH opacities are qualitatively consistent 
with the observed molecular depletion and the chemical 
timescale for L1498; however, the coagulation simulation 
of OH5 and OH8 assumes a density of njj 2 = 10 6 cm' 3 
sustained for 10 5 years (see Ossenkopf 1993). This density 
is substantially higher than the density derived from dust 
continuum models; therefore, the OH5 and OH8 opacities 
may not be based on an accurate representation of fluffy 
aggregation of grains in low density PPCs. 

How do you create a dynamically young yet chemically 
evolved core? If L1498 is in an unstable equilibrium, it 
cannot have been in that state for a long period of time 
since the free-fall time for a 1 — 3 x 10 4 cm" 3 BE sphere, 
tff — yj'in /8Gfimnn (Spitzer 1978), is approximately 
3 — 1.7 x 10 5 years. Clearly, L1498 has not yet collapsed 
and displays chemical differentiation that requires a simi- 
lar timescale. L1498 must have been static or slowly col- 
lapsing for more than 10 5 years. While magnetic fields 
may also play a significant role in cloud support and evo- 
lution for dynamically nascent PPCs (e.g. Mouschovias & 
Spitzer 1976, Li et al. 2002), in the case of L1498, they 
are not necessary for support of the cloud (since Mb e < 
or « Mj). Nevertheless, we can estimate the strength of 
the B-field required for equipartition of magnetic, gravita- 
tional, and kinetic energy. The strength of this B-ficld is 
given by B w 0.51n 1 / 2 Av„ t (km/s) = 8.7 ^G (Lada et al. 
2004). This is a small magnetic field compared with Zee- 
man measurements of the CCS Jn = 3 2 — > 2i transition 
which provide a rough estimate of the linc-of-sight field 
strength of 48 ± 31 [iG (Levin et al. 2001). If L1498 be- 
comes massive enough to be unstable to gravitational col- 
lapse, then the magnetic field appears to be strong enough 
to support L1498, in which case, the timescale for collapse 
is controlled by the ambi-polar diffusion timescale until 
the core becomes magnetically super-critical. 

Recent molecular line studies of PPCs in the Taurus- 
Auriga molecular cloud have identified cores (e.g. L1521B 
and L1521E) that appear to be dynamically young and 
chemically young (Hirota et al. 2002; Tafalla & Santiago 
2004; Hirota, Maezawa, & Yamamoto 2004). The signa- 
ture of chemical youth is strongly peaked sulfur-bearing 
carbon chain molecules (e.g., CCS and C3S), weak NH 3 
and N 2 H + emission, and low molecular depletion of C 18 0. 
Confirmation of the chemical and dynamical youth of these 
cores requires a detailed comparison between submillime- 
ter dust continuum emission and high-resolution observa- 
tions of CCS or C3S. These objects, along with L1498, 
L1512 (n c = 10 5 cm- 3 ), and L1544 (n c = 10 6 cuT 3 ) 
may represent an evolutionary sequence in Taurus- Auriga 
where chemical maturity is reached first and then dynam- 
ical evolution occurs (see Figure 9). 



An important caveat is that our ID modeling has ig- 
nored the observed non-axisymmetry in L1498. The dust 
continuum emission peaks along the north-east ridge and 
not at the FWHM centroid. This effect is most pro- 
nounced in the 1.2mm but is also detected at 850 /im. 
Other PPCs (e.g., L1544 and L63) also display clear non- 
axisymmetries. The spherical BE model is limited and un- 
able to describe this structure. A full 3D radiative trans- 
fer modeling is needed to probe the structure of the non- 
axisymmetry (e.g. Doty et al. 2004). Also, since the 90 
/im ISO image (Ward-Thompson et al. 2002) shows a clear 
gradient across the image, there is evidence that the ISRF 
asymmetrically heats L1498 (also see Langer & Willacy 
2001). 3D models could explore triaxial, BE-like model 
with a spatially varying ISRF compared to fully non- 
axisymmetric density distributions (see Gongalves, Galli, 
& Walmsley 2004) . Encouragingly, a comparison between 
3D models and ID models of L1544 indicate that the ID 
models are an acceptable fit to the average physical struc- 
ture (Doty et al. 2004) 

5. CONCLUSION 

We have presented deep SCUBA observations, BEARS 
and CSO N 2 H + observations, and GBT C3S observations 
of the Pre-protostellar core, Lynds 1498. Radiative trans- 
fer modeling of the submillimeter intensity profiles and the 
complete SED were performed. Our main conclusions are 
as follows. 

L1498 is characterized by a low density Bonnor-Ebert 
sphere with central density of 1 — 3 x 10 4 cm' 3 and outer 
radius greater than 40, 000 AU. The SED is best-fitted 
with Ossenkopf & Hcnning opacities for coagulated grains 
with thick and thin ice mantles (OH8 and OH5 opaci- 
ties). A more realistic treatment of the ISRF was used 
that includes variations in the strength of the incident 
specific intensity and extinction due to dust grains at the 
outer radius of the core. L1498 is a potentially stable 
(Mbe < Mj), magnetically sub-critical core that appears 
to be dynamically young. These results do not change if 
the 1.2 mm continuum peak is used for the radial intensity 
profiles instead of the SCUBA continuum centroid. 

Observations of the depletion of species such as C3S and 
H 2 CO indicate L1498 is chemically evolved (t > 10 5 years). 
The modest N 2 H + abundance is also consistent with a 
chemically evolved source. The best-fitted Ossenkopf & 
Hcnning opacities are consistent with this timescale, al- 
though the Ossenkopf coagulation models assume a higher 
density than is observed toward L1498. The N 2 H + and 
C3S linewidth do not vary significantly with radius out to 
15,000 AU. The N 2 H + abundance also does not vary sig- 
nificantly with radius, while the C3S abundance indicates 
significant depletion toward the center of L1498. 

Comparisons of L1498 with nearby PPCs that have 
also been observed at far-infrared and submillimeter wave- 
lengths indicate that L1498 is larger than average, less lu- 
minous than average, and less massive than average. The 
standard evolutionary indicators used for low-mass proto- 
stars, Tboi and r Lb ° l , do not correlate for PPCs. All of 
the PPCs in this sample are characterized by Tboi < 20 K 
and < 25. 

The temperature and density structure derived from ID 
models may now be used in more realistic models of molec- 
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ular line emission. Future studies of PPCs require high res- such as HAWC on SOFIA and the launch of HERSCHEL 
olution far-infrared observations coupled with 3D radiative will provide the needed high resolution far-infrared obser- 
transfer modeling of multiple wavelengths. Instruments vations. 



APPENDIX 

HYPERFINE LEVELS OF N 2 H+ j = 1 — » 

Since every nucleus in the N 2 H + molecules has non-zero nuclear spin, N 2 H + has a very rich hyperfine spectrum that 
can be exploited to probe the physical conditions in star-forming regions. In this appendix, we summarize the hyperfine 
splitting of the ground state rotational transition of N 2 H+. 

The hyperfine splitting (in order of importance) is due to: (1) an electric quadrupole interaction due to the outer 
nitrogen; (2) an electric quadrupole interaction due to the inner nitrogen; (3) a magnetic dipole interaction due to the 
outer nitrogen; (4) a magnetic dipole interaction due to the inner nitrogen; (5) a magnetic dipole interaction due to the 
hydrogen nucleus; and (6) spin-spin interactions between the nuclei. The final two effects are ignored in calculations of 
hyperfine splitting for astronomical sources since the observed linewidth is > 0.1 km/s (> 31 kHz) and the hyperfine 
coupling constants are small (see Caselli, Myers, & Thaddeus 1995). 

There are 15 electric dipole allowed transitions for N 2 H + J = 1 — > in the \JF\F) basis (AJ = —1, AF\ = 0, ±1, 
AF — 0, ±1, ■/* 0). The 15 transitions are listed in Table 5; however, since in practice the J = level is not split, 
degenerate transitions are collected together into the 7 observed transitions. We adopt the standard notation from Caselli, 
Myers, & Thaddeus for the transition labels (1995). 

It is convenient to calculate the ratio of relative strengths, Ri, of the hyperfine levels for calculations of optical depth 
and column density (§4.2). The relative strengths, s», are determined using irreducible tensor methods (Gordy & Cook 
1984, Chapter 15) and are defined such that the sum of the relative strength of the electric dipole allowed transitions are 
equal to one (see Rudolph 1968). The ratio of relative strengths for the N 2 H+ J = 1 — > transitions in the \JFiF) basis 
are given by 



uifif^ of iF) = 3 -{ \ n M 2 { j f n j 2 n m + 1)i 



(Al) 



allowed v J K J F i =F!,F 1 ,F',F 

F'^F 



where the 6-j symbols needed in Equation Al may be found in Table 5 of Edmonds (1974). The ratios must be summed 
over the degenerate transitions from the same J = 1 F[ F' level (Table 5). 

It is necessary to calculate the dependence of Ri on the linewidth for closely spaced hyperfine lines (see Turner 2001). 
We have performed this calculation for the N 2 H + J = 1 — > transitions assuming each hyperfine component is Gaussian, 



Ej=i *j exp 



-41n2ta^ 



Ri(*v) = — L - ^ . (A2) 

Ej=i^exp(-41n2^) 

The results are shown in Figure 10. For sources with Av < 0.3 km/s, the theoretical hyperfine ratios are sufficient. This 
is the case for L1498. However, for N 2 H + observations of sources with broader linewidths, the lincwidth-corrected Ri 
should be used. 
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Fig. 1. — The jiggle map pointing centers are plotted with the chop direction (a). The (0,0) position is J2000.0 a = 4 h 10 m 52.5 s , 6 = 
+25° 9' 55". Contour maps of L1498 at 170 (b; Ward-Thompson ct al. 2002), 450 fira (c), 850 /im (d), 1200 /im (e; Tafalla et al. 2004), 
and I(T^) of N2H+ J = 1 — > (f). The contour levels are as follows (lowest contour and contour increment in percentage of the peak flux): 
(N 2 H+) 17%(3ct), (1200^m) 20%(3<r) increasing by 20% ; (850Atm) 20%(3<r) increasing by 20% ; (450^m) 25%(3cr) increasing by 25% ; and 
(170^m) 2.0 to 3.6 Jy/pixel increasing by 0.2 Jy/pixel. The 170 image is plotted on a scale that is twice as large as the other 
images. 
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Fig. 2. — Opacity at 350, 450 and 850 fim for August 29 and 30, 1998. The opacity was very stable throughout both nights except for 450 
fim observations during the final two hours of the morning shift (after sunrise). 
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Fig. 3. — The specific intensity of the ISRF with variations in s i3r f and A v (top panel). The dust mass opacities (total gas + dust mass) 
used in dust models (bottom panel). A Mgas/Afdust = 100 was assumed. 
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Fig. 4. — Color coded \r are plotted for the radiative transfer model grid for OH5 dust. x 2 (^85o) an d X 2 (^45o) are plotted for a grid of n c 
and R (top row). x 2 (5i?D) (middle row) and x 2 {^S5o) (bottom row) are plotted for grids of Sj ar f and Ay with outer radii of 20,000, 40,000, 
and 60,000 AU. Note that red points represent the best Xr- 
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Fig. 5. — Dust temperature profile (panel a), SED (b), Is5o(b) profile (c), and /450(b) profile (d) for the best fitted models. The Uranus 
beam profile is shown in the 850 and 450 fim profile. The sector-average radial profiles at 850 fim are plotted as dashcd-lincs (d). 
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Fig. 6. — The standard evolutionary indicators for low-mass protostars are plotted: Tj, ; vs. L boL . 



PPCs are plotted as circles and Class 



sources (Shirley et al. 2000, Young ct al. 2003) arc plotted as triangles. There is no correlation between these indicators for PPCs; however, 
all observed PPCs have T bol < 20 K and , Lhal < 25. 
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Fig. 7.— The N 2 H+ J = 1 -> spectrum (a) towards the SCUBA continuum centroid (J2000.0 a = A h 10 m 53.3 s , <5 = +25° 9' 39"). 
The lincwidth of the isolated hyperfine component (JF\F = 101 — » 012) is plotted versus radius (b). The N2H + J = 3 — + 2 spectrum at the 
SCUBA continuum centroid is shown in panel (c) with a theoretical hyperfine stick spectrum. The N2H+ J = 3 — > 2 detection (Crapsi et al. 
2004) toward the 1.2 mm continuum peak (J2000.0 a = 4 h 10 m 51.5 s , <5 = +25° 9' 58") is shown in the inset (c). The abundance of N 2 H+ 
versus radius is plotted in panel (d) with statistical errorbars (07, cr T , & cs^)- No significant variations of the linewidth and abundance are 
observed on scales out to 15000 AU. 
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Fig. 8. — The C3S J = 4 — > 3 integrated intensity map is plotted in panel (a) with a cross at the position of the SCUBA dust continuum 
ccntroid. The (0,0) position is identical to Figure 1 (J2000.0 a = 4 h 10 m 52.5 s , <5 = +25° 9' 55")- The variation of linewidth is shown with 
radius (b). The C 3 S and C 3 34 S (inset) J = 4 -> 3 spectra (c) towards the SE C 3 S peak (J2000.0 a = A h 10™ 56.9 s , <5 = +25° 9' 22'.'5 ). The 
abundance of C3S versus radius in the major axis sector (§2.1.2) is plotted in panel (d) with statistical errorbars (07 & cs„)- ^? = is the 
dust continuum ccntroid position; — R refers to positions to the SE. Significant depletion is seen toward the dust continuum ccntroid. 
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Evolutionary Sequence 9 

Fig. 9. — Is this a possible evolutionary sequence for PPCs in Taurus-Auriga? The 850 fim images of L1512 and L1544 are from Shirley et 
al. (2000) while the modeled central densities are from Evans et al. (2001). Lee et al. (2003) have characterized the chemical states of L1512 
and L1544. There arc no published or archive submillimeter detections of the nascent pre-protostellar candidate, L1521B. 
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Fig. 10. — The variation of the ratio of relative strengths for N2H + J = 1 — ► hypcrfine lines with lincwidth. The theoretical ratios are 
valid for Av < 0.3 km/s. 
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L1498 Spectral Energy Distribution 
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References. - 1. Ward- Thompson 
et al. 2002, 2. Ward-Thompson et al. 
1994, 3. This paper, 4. Tafalla et al. 
2004 



Table A2 



Dust Opacity Properties 



Model 
Name 


Description 


a a 
Ksmm 


k„(850) 
(cm 2 g^ 1 ) 


Ref. 


Draine-Lee 


Silicates + Carbonaceous grains 


2.0 


0.003 


1 


WD3 


Silicate + Carbonaceous + PAH ISM fit 


1.8 


0.005 


2 


WD4 


Silicate + Carbonaceous + PAH ISM fit 


1.6 


0.006 


2 


Pollackl 


Silicates + Organic C 


2.2 


0.004 


3 


MMP 


Empirical Fit to ISM 


2.3 


0.030 


4 


OH2 


Coagulated for 10 5 yrs. no ice mantles 


1.3 


0.034 
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OH5 


Coagulated for 10 5 yrs. thin ice mantles 
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OH8 


Coagulated for 10 5 yrs. thick ice mantles 


1.9 


0.022 
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a K <~ where /3 smm is determined by a linear regression over A € [350/xm, 1300^m]. 



References. — 1. Draine & Lee 1984, 2. Weingartner & Draine 2001, 3. Pollack et 
al. 1994, 4. Mathis, Metzger, & Panagia 1983, 5. Ossenkopf & Henning 1994 



Table A3 
Best-fitted Models 3, 
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b Model is consistent with Cambresy 1999 optical extinction constraint that 
Ay < 3 mag. 



Table A4 

Properties of Nearby PPCs Observed with ISOPHOT a AND SCUBA 
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0.11 


7.4 
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1830 


0.4 


0.09 


6.0 
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a ISOPHOT fluxes from Ward-Thompson et al. 2002. 

b SCUBA 850 (im continuum FWHM contour centroid in J2000.0 coordinates. 
c Reference for SCUBA observations. 

dm _ hC(4) l^vSudv 
J-bol — 4fcc(5 ) joo Svdv ■ 

c Dust-determined mass in a 120" aperture. 
f Median of the sample. 



References. — 1. This paper; 2. Shirley et al. 2000; 3. CADC SCUBA Archive; 4. Bianchi et al. 2003, Lai et al. 2003 
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Table A5 

N 2 H + J = 1 — > Hyperfine Transitions 

J'F[F' -» JFiF J'i^F' -» JFi^ 1 v Au i?, 

Standard Notation (GHz) (km/s) 

1 1 -> 1 1 ^ 110^011 93.1716200 +6.944 | 

112^012 | 112^ 1 2 93.1719168 +5.988 § 

m-012 $ 

m-011 3 

111^010 ^ 111^ 1 93.1720533 +5.549 | 

12 2-^012 J 

1 2 2 -> 1 1 ^ 122^011 93.1734796 +0.956 § 

123^012 g 123^012 93.1737767 +0.000 1 

121^012 ^ 

121->011 j|g 121-> 1 1 93.1739666 -0.611 f 

1 2 1 -> 1 J, 

101^012 J 101^ 1 2 93.1762650 -8.011 f~ 

101-tOll I 

1 1 -> 1 %r 



